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ABSTRACT 


The  problem  of  the  flow  of  a gas-liquid  mixture  through  a multi- 
stage axial  compressor  originally  designed  for  air  flow  arises  during 
take-off  from  a rough  runway  with  water  on  it  and  during  rain.  The 
problem  is  complicated  because  of  the  presence  of  discrete  droplets  in  a 
rotating  compressor  flow  field  and  because  of  the  phase  change  during 
compression.  In  order  to  alleviate  this  problem,  it  is  necessary  to 
examine  both  changes  in  blade  and  casing  design  as  well  as  possible 

additional  controls  including  variable  geometry.  The  objectives  of  the 

/ 

investigation  presented  here  are  (a)  to  establish  estimates  of  possible 
changes  in  performance  of  typical  aircraft  compressors  operating  with 
two-phase  flow,  (b)  to  determine  the  governing  aerothermodynamic  equations 
for  two-phase  flow  in  compressors,  and  (c)  to  outline  directions  in  which 
further  research  iray  prove  fruitful  in  this  subject.  Preliminary 
investigations  have  revealed  the  problem  areas  in  the  fan,  L.P.,  and  H.P. 
compressor  stages.  The  basic  aerothermodynamic  equations  have  been 
deduced  in  a form  suitable  for  considering  development  of  a computational 
program. 
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1 . INTRODUCTION 


Two-phase  flow  in  turbomachinery  is  of  interest  in  a number  of 
technological  problems,  e.g.  rain  ingestion  into  compressors  on  ground  or 
during  flight,  condensation  in  steam  turbines,  pumping  of  crude  oil, 
circulation  of  coolants  in  nuclear  reactors,  etc.  The  principal  interest 
in  this  Report  is  in  the  problem  of  air-water  mixture  flow  through 
diffusers  and  compressors.  In  gas  turbines  installed  in  aircraft,  the 
problem  of  water  ingestion  into  the  engine  is  known  to  lead  to  serious 
consequences,  both  by  reducing  the  stall  margin  and  occasionally  leading 
to  after-burner  blowout  (Ref.  1).  The  ingestion  of  water  may  arise 
during  take-off  from  puddles  on  rough  runv/ays  or  during  flight  through  a 
rainfall.  It  is  unclear  at  the  moment  if  the  problems  of  surge  and  blow- 
out can  be  removed  with  simple  bleed-type  controls  or  if  one  has  to 
employ,  for  instance,  casing  treatment  or  even  redesign  of  selected 
blading  in  extreme  cases.  It  is  possible  that  no  general  solution  can  be 
obtained  for  all  types  of  compressors  and  all  types  of  water  ingestion, 
but  one  may  be  able  to  establish  general  guidelines  for  design,  control, 
and  operation  of  compressors.  Such  general  guidelines  may  then  be 
adapted  to  suit  given  compressor  configurations. 

No  adequate  theory  exists  for  the  motion  of  a two-phase  fluid  through 
a compressor.  The  compressor  of  interest  in  this  Report  is  a rotating 
machine,  for  e.g.,  an  axial  flow  compressor.  The  increase  in  pressure 
arises  on  account  of  work  done  on  the  fluid.  In  the  absence  of  rotation 
and  work  input  to  the  fluid,  one  can  consider  a simple  diffuser  with  an 
adverse  pressure  gradient.  There  is  no  adequate  theory  for  two  phase 
flow  even  through  such  a diffuser 
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The  two-phase  flow  of  interest  in  this  Report  has  two  basic 


complexities: 

(a)  In  the  air-water  mixture  flow,  water  is  in  the  form  of  discrete 
droplets  which  may  undergo  break-up  or  coalescence  during 
motion. 

(b)  When  work  is  done  on  the  fluid,  water  may  undergo  a phase  change, 
leading  ultimately  to  the  fluid  becoming  an  air-steam  mixture. 

The  aerothermodynamics  of  two-phase  flow  thus  becomes  quite  complicated. 

Considering  an  axial  compressor,  such  as  for  e.g.  the  TF-30,  several 
additional  complexities  arise  in  the  flow: 

(a)  Alternating  rotating  and  stationary  blade  rows. 

(b)  Division  of  the  compressor  into  the  L.P.  and  H.P.  sections. 

(c)  Flow  entry  geometry  at  the  inlet  to  the  compressor. 

Taken  together  with  the  complexities  of  two-phase  flow,  the  foregoing 
presents  a formidable  problem. 

At  this  time,  there  are  no  experimental  studies  of  a sufficiently 
detailed  nature  to  provide  some  guidance  for  theory  and  analysis.  The 
experimental  studies  again  prasen*  complicated  problems  on 

account  of  the  rotation  in  the  fluid  (both  centrifugal  and  Coriolis 
forces  have  to  be  included  appropriately)  and  the  presence  of  the  dis- 
crete phase.  For  example,  the  combined  effects  of  centrifuging  and 
coalescence  can  lead  to  a film  flow  of  water  along  the  casing  of  a 
compressor  and  it  is  important  to  establish  the  possibility  of  such  a 
film  flow  both  from  an  analytical  and  a practical  point  of  view. 
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An  interesting  aspect  of  the  presence  of  water  at  entry  to  the 
compressor  is  the  - -ase  in  humidity  of  air  throughout  the  compressor. 


r 


A measure  lo-  ' •v.'yif-y  is  what  may  be  termed  the  "Saturation  Mass 
Fraction  of  Wate:  , - .fined  as  the  mass  of  water  in  a pound  mass  of  air 
at  saturation  under  given  conditions  of  pressure  and  temperature.  Thus, 
if  Wg  is  the  satur  ^ mass  fraction  of  water  at  given  values  of  ambient 
temperature  and  pre.  sure,  one  can  write  (Ref.  2) 

Ws  = 0.6219  PS/(P  - Ps)  0.1) 

where 

Molecular  Height  of  Hater  Vapor 
~ Molecular  Weight  of  DrJ  Air 

P = air- water  vapor  mixture  pressure 

= partial  vapor  pressure  of  water  at  conditions  of 
saturation. 

W$,  by  definition,  is  a local  parameter.  If  one  considers  standard 
atmospheric  conditions  at  entry  to  the  P&W  TF-30  compressor  ano  the  local 
conditions  at  the  exit  of  stages  1,  2,  3,  4,  5,  and  6,  one  can  tabulate 
the  values  of  W£  at  the  local  conditions  as  shown  in  Table  1.1.  Thus, 
at  a pressure  of  39.71  psia  and  temperature  of  419.78°K,  the  value  of  W 
is  0.9821,  which  means  that  the  amount  of  water  that  can  be  hold  in  1 lb 
of  air  under  saturated  conditions  is  0.9821  lb. 

One  snould  hov/ever  not  be  misled  by  the  possible  values  of  W . In 
actual  practice,  the  attainment  of  equilibrium  conditions  with  respect  to 
humidity  requires  a finite  length  of  time  which  in  turn  depends  upon  the 
size,  number  and  geographical  distributions  of  droplets,  and  their  convec- 
tion velocity.  There  is  no  clear  method  of  establishing  the  increase  in 
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Table  1 .1 


Local  Values  Of  W Under  Saturated  Conditions* 
s 


*Fan  Data  Taken  At  Diameter  Equal  to  Outer  Diameter  Of  L.P.C.  Capture 
Area.  L.P.C.  Data  Taken  At  Blade  Tips. 
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humidity  of  the  air  in  the  presence  of  water  droplets  over  a unit  period 
of  time  even  in  a simple  duct  flow.  In  a complicated  flov/  such  as  that 
in  a compressor,  the  problem  becomes  practically  insurmountable.  However, 
one  has  to  allow  for  some  change  in  the  liquid  water  content  along  the 
compressor  on  account  of  changes  in  the  humidity  of  air. 

The  objective  of  this  Report  is  three-fold:  (1)  To  conduct  a 

performance  estimation  study  on  a selected  compressor  for  the  general 
effects  of  two-phase  flow  in  the  compressor  when  it  has  been  originally 
designed  for  air  flc<v  only;  (2)  To  establish  the  basic  governing  equations 
for  the  aerothermodynamics  of  tv/o-phase  flow  through  a compressor;  and 
(3)  To  outline  a program  of  theoretical  and  experimental  studies  that  can 
lead  to  an  engineering  solution  to  the  problem  of  water  ingestion  into 
installed  engine  compressors. 


1.1  Outline  of  Report 

In  order  to  set  the  background  for  the  problem  of  two-phase  flow 
through  a compressor  we  present  in  Section  2 some  basic  two-phase  flow 
considerations.  The  presence  of  water  in  discrete  phase  introduces  in  a 
compressor  problems  of  heat  and  mass  transfer  and  of  phase  change.  We 
therefore  discuss  thermo-physical  and  transport  considerations. 

In  Section  3 some  preliminary  results  are  presented  on  the  changes 
in  the  overall  performance  of  selected  compressors  when  two-phase  flow  is 
introduced  into  them.  A general  discussion  on  aerothermodynamics  of  two- 
phase  flow  through  a compressor  is  presented  in  Section  4. 

Finally,  an  outline  for  future  studies  is  developed  in  Section  5 
with  a special  emphasis  on  experimental  studies. 
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2.  TWO-PHASE  FLOW  CONSIDERATIONS 
The  subject  of  two-phase  flow  is  of  interest  in  many  technologies 
and  several  reviews  of  the  subject  are  available  (Refs.  3-7)  although  the 
subject  is  extremely  complicated  and  presents  many  unsolved  problems.  We 
shall  concentrate  here  on  those  aspects  of  two-phase  flow  that  are 
specialized  to  the  current  investigation. 

The  two-phase  flow  of  interest  in  this  investigation  is  air  carrying 
water  droplets.  The  water  droplets  are  the  result  of  rain  ingestion  or 
of  water  sprayed  into  the  inlet  during  take-off  from  puddles  on  uneven 
runways . 

The  presence  of  water  alters  the  humidity  of  air  at  entry  and  in 
fact  throughout  the  compressor  with  changes  in  temperature  and  pressure 
so  long  as  water  is  present  in  the  flow 

The  first  consideration  when  water  is  present  in  discrete  droplet 
form  is  the  size  and  number  density  distribution.  If  the  water  is 
ingested  on  account  of  rain,  the  initial  size  and  number  density  distribu- 
tion can  be  related  in  a very  approximate  fashion  to  the  rain  droplet 
size  distribution  (Refs.  8-9).  When  water  is  splashed  into  the  compressor 
inlet  from  runway  puddles,  one  can  make  a preliminary  estimate  of  the 
amount  of  water  that  may  be  splashed  by  the  wheel  on  account  of  the 
relative  motion  between  the  water  surface  and  the  adjoining  air  flow 
(see  Appendix  I).  There  is  considerably  greater  ambiguity  in  the  spray 
I characteristics  (Ref.  10).  In  addition  to  the  uncertainties  in  the  spray 

I 

envelope  from  the  wheel,  there  is  also  the  difficulty  of  determining  the 
influence  of  inlet  geometry  and  flow  on  the  actual  ingestion  process. 
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The  inlet  flow  itself  is  altered  in  many  cases  by  the  proximity  of  the 
ground.  Until  such  time  when  better  data  become  available  on  droplet 
size  and  number  characteristics,  the  only  feasible  approach  is  to  assume 
that  distributions  of  statistical  mean  diameter  and  number  density  may  be 
treated  as  parameters  at  the  inlet  to  the  compressor. 

Upon  entry  ihto  the  compressor,  the  size  and  number  density 
distributions  of  droplets  alter  continuously  on  account  of  the  foil  ow  ing. 
(i)  break-up  and  coagulation;  and 
( i i ) evaporation. 

The  break-up  and  coagulation  arise  on  account  of  the  motion  in  a moving 
and  generally  rotating  environment  in  the  presence  of  gravitational 
action.  The  evaporation  requires  the  droplets  to  become  heated  to  the 
boiling  point  at  the  local  pressure  and  then  to  evaporate  with  the 
absorption  of  latent  heat.  The  evaporation  is  a strongly  rate  dependent 
process  and  several  characteristic  length  and  time  scales  are  involved. 

At  this  stage  the  only  feasible  approach  to  take  in  such  problems  is  to 
include  the  break-up  and  evaporation  processes  in  the  changes  in  the 
statistical  mean  diameter  and  number  density  locally  along  the  streamlines. 

The  second  consideration  when  water  is  present  in  discrete  droplet 
form  is  the  distribution  of  velocities  for  the  different  droplets.  We 
assign  a single  mean  velocity  for  each  group  of  droplets  with  an  assigned 
mean  statistical  diameter.  We  are  thus  avoiding  consideration  of  the 
motion  of  individual  droplets  in  any  group. 

The  third  consideration  arises  from  the  fact  that  water  droplets 
introduce  the  drag  of  droplets  into  the  momentum  balance  of  the  air  flow. 
The  drag  is  based  here,  by  assumption,  on  the  Reynolds  number  calculated 
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with  respect  to  the  statistical  mean  diameter  and  the  relative  velocity 
of  the  droplets  with  respect  to  air  flow  velocity.  The  drag  of  N$ 
droplets  is  taken  as  Ns  times  the  drag  of  one  dg  size  particle. 

The  fourth  and  last  consideration  to  be  taken  into  account  here  is 
the  difference  between  the  compressibility  of  the  mixture  and  the 
compressibility  of  air.  The  acoustic  velocity  in  the  mixture  is 
different  from  that  in  the  air  and  the  local  Mach  number  of  the  mixture 
at  any  location  is  therefore  based  on  the  acoustic  velocity  in  the 
mixture. 

In  view  of  the  foregoing  considerations,  some  tentative  bases  have 
been  selected  to  take  into  account  the  various  processes.  They  are 
elucidated  in  the  following. 

2.1  Water-Phase  Continuum 

The  water  content  in  the  mixture  is  treated  as  a continuum  in  the 
general  equations  of  motion  along  with  air  and  the  mixture.  The  water 
content  in  a unit  volume  at  any  location  is  assigned  a value  oy  by 
volume.  Continuity  is  thus  assumed  for  (1  - oy)  of  air  and  <yy  of 
water.  It  should  be  pointed  out  that  changes  along  the  flow  and 

hence  is  a dependent  variable  in  the  flow  equations. 

In  spite  of  the  assumption  of  continuum  for  the  water  droplets,  it 
is  assumed  that  the  water  droplets  do  not  contribute  to  the  local 
pressure  based  on  the  reasoning  that  the  droplets  are  not  distributed 
randomly  enough. 

It  may  be  pointed  out  here  that  for  certain  purposes,  the  droplets 
are  considered  as  discrete  entities.  Three  processes  for  which  it  is 
considered  necessary  to  discretize  the  droplets  are  (1)  the  exchange  of 


momentum  between  the  air  and  the  droplets,  (2)  the  evaporation  of  the 
droplets,  and  (3)  the  break-up  and  agglomeration  of  droplets. 


2.2  Compressibility  of  the  Mixture 

The  acoustic  velocity  in  a droplet-laden  air  flow  is  given  by  the 


following  relation  (Ref.  11). 


c = jO  - <\>g  + av  p 


1 - a 


V c 2 


(2.1) 


g g cp  cp  J 


where 


c = mixture  sonic  velocity 

c?v  = particulate  liquid  volume  fraction  = volume  of  particulate 
liquid  per  unit  volume  of  mixture 
p = densi ty 

( ) = pertaining  to  gas  phase 
( )p  = pertaining  to  particulate  liquid  phase. 

Equation  (2.1)  is  presented  graphically  in  Figure  2.1  for  standard 
atmospheric  conditions  and  for  a range  of  cy  from  0.00  to  0.10;  the 
range  over  which  this  report  is  interested. 


2.3  Distribution  of  Droplet  Size  and  Humber  Density 

The  volume  oy  of  water  in  a unit  volume  of  the  mixture  at  any 
location  is  specified  by  the  following  relation. 


; = r dj  H 
v 6 s s 


(2.2) 


where 


= statistical  mean  droplet  diameter  of  a cloud  of  liquid  droplets 
- statistical  droplet  number  density  of  a cloud  of  liquid 


droplets  with  a statistical  mean  diameter,  d 
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For  our  purposes  the  statistical  mean  droplet  diameter  is  taken  as 
the  Sauter  mean  diameter  (Ref.  12) 

d,  ■ (d"aX  d3p(d)dd/fd"’aX  d2p(d)dd  (2.3) 

Jo  Jo 

where 

d = an  individual  droplet  diameter 

p(d)  = the  probability  of  a droplet  having  a diameter  between  d 
and  (d  + <$d). 

It  should  be  noted  that  a droplet  with  diameter  equal  to  the  Sauter  mean 
diameter  has  the  same  surface  area  to  volume  ratio  as  the  entire  spray. 

It  should  be  recognized  that  in  addition  to  oy,  a$  and  N$  also  may 
change  along  the  flow  depending  upon  droplet  break-up  and  evaporation 
processes.  Thus,  when  av  changes,  there  may  also  arise  a change  in  d^ 
and  then  one  can  assign  a value  for  Ns  corresponding  to  the  changed 
values  of  oy  and  ds. 

At  entry  to  the  compressor,  a radial  and  circumferential  distribu- 
tion of  oy  and  d$  is  specified,  although  under  an  axisymmetric  assumption 
the  circumferential  variation  in  distribution  can  be  neglected. 

2.4  Break-up  of  Droplets 

A droplet  is  assumed  to  breakup  when  the  Weber  number  becomes  of  the 
order  of  ten.  The  Weber  number  based  on  dg  is  defined  by  the  relation 
(Ref.  13): 

We  = Pg(Vg  - Vp)2ds/a  (2.4) 
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where 

We  = Weber  number 
V = velocity 

a - droplet  surface  tension. 

It  may  be  noted  again  that  is  not  affected  by  the  break-up 
locally,  although  the  change  in  ay  from  position  to  position  is  affected 
by  the  droplet  break-up. 


2 . 5 Formation  of  a Film  of  Water 

The  combined  effect  of  droplet  agglomeration,  centrifugal  action  and 
gravitational  action  can  lead  to  the  formation  of  a stream  of  water  at 
the  casing  wall.  There  is  no  experimental  evidence  for  this  and  it  is  a 
very  important  question  both  from  the  point  of  view  of  fluid  dynamics 
and  the  possible  arrangements  required  for  the  removal  of  water.  It  is 
extremely  difficult  tc  establish  the  formation  of  a film  of  water  on  an 
analytical  basis. 

We  assume  that  when  the  value  of  c?y  reaches  a large  enough  value 
for  the  mixture  to  be  treated  as  a bubbly  mixture,  there  arises  a film. 
Such  a value  of  is  ustimated  to  be  0.40  for  tap  water  and  air  and 
0.90  for  "pure"  water  and  air  (Ref.  14)  at  a pressure  of  one  atmosphere. 


2.6  Drag  on  a Cloud  of  Droplets 

The  drag  on  a spherical  droplet  carried  in  an  airflow  is  a function 
of  its  Reynolds  number  based  on  the  droplet  diameter  (dg)  and  the  droplet 
velocity  relative  tG  that  of  air  velocity  in  the  vicinity  of  the  droplet. 
If  we  assume  de  = 1.0  mm  and  a relative  velocity  of  10  m/sec.,  the 
Reynolds  number  becomes  688  at  standard  atmospheric  temperature  and 
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pressure.  The  drag  on  a single  droplet  may  then  be  written  as  follows 
(Ref.  15). 

° = °-5  w v v\  (2-5) 

where 


CD  = 24.0  Rep-‘-° 

CD  = 27.0  Rep'°‘84 
CQ  = 0.414  Rep0-11433 


Rep  < 0.5 

0.5  < Re„  < 70.0 
P 

70.0  < Rep  < 1.3  x 106 
1.3  x 106  < Rep. 


The  particulate  Reynolds  number  is  defined  as  follows. 


(2.6) 


where 


d 

s 


- V. 


(2.7) 


CD  = particle  drag  coefficient 

Ap  = projected  frontal  area  of  an  equivalent  sphere  having  the 
same  mass  as  the  droplet. 

The  drag  on  a cloud  of  droplets  becomes  much  more  complicated  to 
estimate  because  of  (a)  the  mutual  interaction  between  the  droplets  and 
(b)  the  arbitrariness  of  the  direction  of  motion  of  different  droplets 
at  different  times  in  the  vicinity  of  any  one  droplet.  References  16-28 
provide  a short  bibliography  on  this  subject.  Two  examples  of  drag 
correlations  are  given  here,  one  due  to  Refs.  16  and  21,  and  the  other 
due  to  Ref.  22. 


2.6.1  Sop's  Drag  Coefficient  Using  Ergun's  Pressure  Drop  Equation 

Drag  of  a cloud  of  N droplets  = N_  x the  drag  of  one  dr  size  particle 

^ b s 

* Ns  x °'5  CD  pg<vg  - yV  <2'8> 
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where 


A = projected  frontal  area  of  a droplet  having  the  statistical 
ps 

mean  diameter  d . 


200  cv  |jq  7 

- (1  - cv)ds  u'  PS  3(1  - ov] 


(2.9) 


0.08  < a.,  < 0.60  and  1.0  < — < 3000.0 
V °v 

u - absolute  viscosity 

uE  - relative  superficial  velocity/(l  - av). 

The  relative  superficial  velocity  is  the  relative  gas  flow,  (V  - V ), 

y r 

based  upon  the  unobstructed  flow  area,  i.e.,  the  flow  area  when  no 
droplets  are  present. 


2.6.2  Tam’s  Formula  for  the  Drag  Exerted  on  a Particle  in  a Cloud  of 


Spherical  Particles 

Drag  on  one  d$  size  particle  = D(a)  U0 


(2.10) 


where 


U0  e = the  mean  fluid  velocity  before  the  addition  of  the  d size 

A j 

particle 

D(o)  = a drag  coefficient  which  depends  on  the  size  of  the  d$ 

particle  and  the  statistical  properties  of  the  N$  particles 


D(a)  = 3 it  d$[l  + ~2~  + ds)3 

3a  % 

6t?  m2  + [30ir2  m2  + 24ir  m^l  — jp-)] 
a = — 2 - 3a.. 


■ 

m(n)  = n(ar)arn  da, 


(2.11) 


(2.12) 


(2.13) 


and  uihoy-A 

Oil U W 1 1C  I C 

a = particle  radius 
n(or)  = particle  distribution  function 
m(n)  = the  nth  moment  of  the  particle  distribution  function 
0.0  < oy  < 0.50  0.0  < Rep  < 2.0. 

2.7  Heat  (Mass)  Transfer  Between  a Cloud  of  Droplets  and  the  Surrounding 
Fluid 

In  our  case  of  interest,  the  interfacial  transfer  of  heat  and  mass 
betv/een  a single  droplet  and  the  surrounding  fluid  is  primarily  due  to 
forced  convection.  Free  convection  and  radiative  (heat)  transfer  are 
neglected.  The  transfer  rates  may  be  determined  by  the  following. 

Heat  Transfer  Rate  = h^  A(T  - T ) 

(2.14) 

Mass  Transfer  Rate  = hm  A(C,  - C ) 

m ws  wg 

where 

= average  heat  transfer  coefficient  for  the  droplet  (cloud) 

~ average  mass  transfer  coefficient  for  the  droplet  (cloud) 

A - droplet  (cloud)  surface  area 

T = temperature 

Cws  = water  vapor  concentration  at  surface  droplet  (cloud) 

CWg  = water  vapor  concentration  in  fluid  flowing  around  droplet 
(cloud) . 

When  the  relative  velocity  between  a single  droplet  and  the 
surrounding  fluid  approaches  zero  the  following  relationship  is  used  to 
determine  the  transfer  rates  (Ref.  29). 

Nil  = 2.0  ] 
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where 


i 

j 


Nu  = Nusselt  number  = h.  d /k 

n s 

Sh  = Sherwood  number  - h d„/D. 

m s b 

k = thermal  conductivity  of  the  fluid  surrounding  the  droplet 
(cloud) 

= binary  mass  diffusivity. 

The  convective  interfacial  heat  and  mass  transfers  are  a function  of 
the  respective  Peclet  number,  and  the  internal  circulation  of  the  droplet. 
The  Peclet  numbers  are  defined  as  follows: 


Peh  = Re„  x Pr 
h d 


dr/«. 


(9 

\ “ * 1 w / 


Pem  = Rep  x Sc  = 


V - V 

g p 


ds/°b 


where 


Peh  = heat  transfer  Peclet  number 

Pem  = mass  transfer  Peclet  number 

Pr  = Prandtl  number  = ug/pg 

Sc  = Schmidt  number  = y /p  D, 

y y u 


(2.17) 


a = thermal  diffusivity  of  the  fluid  surrounding  the  droplet 
(cloud). 

Again  assuming  d$  = 1 ram  and  a relative  velocity  of  10  m/sec,  Pe^  = 498, 
and  Pe^  = 391  at  standard  atmospheric  conditions. 

When  the  particulate  phase  changes  from  discrete  droplets  to  bubbles 
the  internal  circulation  varies  from  zero  to  a maximum.  Fluid  droplets 
are  considered  to  have  an  intermediate  amount  of  internal  circulation. 

The  presence  of  surfactant  (surface  reactant)  impurities  contaminates 
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the  gas-fluid  particle  interface  and  also  tends  to  reduce  the  internal 
circulation  in  droplets  and  bubbles.  This  in  turn  considerably  reduces 
the  rates  of  interfacial  transfer. 


References  30  and  31  provide  a summary  of  a large  number  of  analytic 
and  correlated  results  for  single  particle  interfacial  transfer.  Three 
such  results  follow. 

Friedlander  (Ref.  32). 

Pe.  Pe.  2 

Nu  = 2[1  + + -ys—  + ...  ] 


■ ■ 4^  • 12 

Pem  Pe  2 

Sh  = 2[1  + —jj~  + + 


(2.18) 


These  results  apply  for  no  internal  circulation,  small  Peclet  numbers, 
and  small  particle  Reynolds  numbers. 

Levich  (Ref.  33). 


Nu  = 0.998  Peh/3  + 2 


Sh  = 0.998  Pe'3  + 2 
m 


(2.19) 


Levich' s results  apply  for  no  internal  circulation,  large  Peclet  numbers, 
and  small  particle  Reynolds  numters. 

Ruckenstein  (Ref.  34). 


where 


Nu  = 0.895 (I)3*  Pe^  + 2 
Sh  = 0.895(|)*  PeJ*  + 2 


= 2(1  + 6) 


<2. 20) 


(2.21) 


6 = ug/(up  + 5) 


(2.22) 


e = "interfacial  viscosity"  due  to  adsorbed  surfactant  impurities. 
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These  results  apply  for  mode rate  to  strong  internal  circulation,  large 
Peclet  numbers,  and  small  particle  Reynolds  numbers. 

When  the  relative  velocity  between  a cloud  of  droplets  and  the 
surrounding  fluid  approaches  zero  the  following  relationship  may  be  used 
to  deternrne  the  transfer  rates  (Ref.  31). 


Nu  = 2/(1 


) 1 


(2.23) 


Sh  = 2/(1  - ovl/3)  J 

The  Mussel t and  Sherwood  numbers  are  defined  in  the  same  manner  as  for 
the  single  droplet  case,  but.  note  that  the  fluid  surrounding  the  cloud  of 
particles  may  well  include  other  particle  clouds  and  therefore  the 
thermal  conductivity,  thermal  diffusivity,  and  binary  diffusion  coeffi- 
cient must  be  modified  accordingly. 

The  interfacial  convective  heat  (or  mass)  transfer  between  a cloud 
of  droplets  and  the  surrounding  fluid  is  not  only  a function  of  the 


Peclet  number  (Pe^ »Pem)  and  the  internal  circulation  of  the  particles. 

It  is  also  a function  of  particle  concentration  and  particle  size 
distribution.  These  reflect  the  influence  of  both  particle-particle 
interaction  and  hydrodynamic  fluid-particle  interaction.  References 
(35-44)  provide  a short  bibliography  on  this  subject.  Two  examples  of 
correlations  are  given  here,  one  due  to  Reference  36  and  the  other  due  to 
Reference  41 . 

Pfeffer  (Ref  36). 


Nu  = 1.26 


Sh  = 1.26 


(2.24) 
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where 

V = 2 + 28  + av5/3(3  - 2g)  (2.25) 

X = 3 + 26  + 2a  ^(1  - 6).  (2.26) 

Note  that  the  relative  velocity  used  to  define  the  cloud  Peclet 
number  is  defined  as  the  difference  between  the  gas  flow  velocity  and  the 
mean  velocity  assigned  to  the  cloud  of  particles.  Refer  to  the  second 
consideration  listed  in  Section  2. 

Pfeifer's  results  apply  to  the  flow  of  an  ensemble  of  uniform  size 
particles  with  no  internal  circulation,  a large  Peclet  number,  and  a 
small  particle  Reynolds  number. 

Yaron  and  Gal -or  (Ref.  41) 


Nu 


Sh 


1.1284  Pe^ 

o - °,>* 


n 


1.1284  Pe* 

o : %)* 


[1 


2.8852 


2.8852 

Rep* 


F(u  ,p  ) ] 
F(y,p)] 


(2.27) 

(2.28) 


where  K(e.p)  = [1  + ^/[l  + 

g g 

These  results  apply  to  the  flow  of  a cloud  of  uniform  size  particles 
with  no  restriction  on  the  amount  of  internal  circulation,  a large  Peclet 
number,  and  50  < Re  < 1500. 

r 

As  previously  mentioned,  we  are  interested  in  the  variation  of  the 
thermal  conductivity  of  the  overall  mixture  with  particulate  volume 
fraction,  o . Gotoh,  Reference  45,  recommends  without  restrictions  the 
use  of  the  following  equation  for  evaluating  the  effective  thermal 
conductivity  of  two-phase  heterogeneous  flow. 
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(2.29) 


N 


where 

k = the  effective  thermal  conductivity  of  the  mixture, 
e 

Yamada  and  Takahashi,  Reference  46,  have  conducted  an  interesting 
study  on  the  variation  of  the  effective  thermal  conductivity  of  suspen- 
sions with  particle  shape,  orientation,  dimension,  and  surface  condition. 
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3.  PERFORMANCE  ESTIMATION 


In  light  of  the  discussion  in  Section  2,  a multi-stage  axial 
comDressor  operating  with  a two-phase  fluid  is  studied  best  by  dividing 
it  into  three  parts  as  follows: 

1)  The  low  pressure  stages  in  which  no  substantial  change  of  phase 
(from  liquid  to  vapor)  has  occured; 

2)  The  intermediate  stages  in  which  there  is  a complete  change  of 
phase;  and 

3)  The  high  pressure  stages  in  which  one  has  essentially  an  air- 
vapor  mixture. 

In  each  group  there  is  a change  of  performance,  efficiency,  and  surge 
margin. 

In  view  of  the  uncertainties  on  transport  properties  of  particulate 
two-phase  flow,  we  have  concentrated  in  the  present  investigation  on  the 
performance  of  a compressor  operating  with  air- vapor  mixture.  This 
effort  is  divided  into  the  following: 

(1)  Performance  calculation  of  a single  stage  air  compressor, 

PURDUE  001  & 002  Compressors 

(2)  Performance  calculation  for  the  multi-stage  compressor,  N.6.T.E. 
$1 09  Compressor 

(3)  Performance  calculation  for  the  P&W  TF-30  Compressor. 

Each  of  the  fo rementioned  are  evaluated  using  temperatures  corresponding 
to:  (1)  design  conditions  - 288  K,  (2)  air  containing  water  vapor  in  the 

form  of  dry  steam  - 373  K,  and  (3)  air  before  the  latent  heat  is  absorbed 
due  to  evaporation  - 429.04  K or  484.79  K. 


V*I 


3.1  Single  Stage  Compressor 

A single  stage  compressor  was  desinged  with  the  following  design 


parameters: 


mass  flow 


20  kg/sec 


stage  temperature  rise  = 2Q°C 


r.p.m. 

mean  blade  speed 


= 9,000 
= 180  m/sec 


Both  a constant  reaction  (PURDUE  001)  and  a free-vortex  (PURDUE  002) 
compressors  were  designed. 

The  performance  of  the  single  stage  compressor  was  calculated  for 
operation  with  the  following  working  fluids: 

a)  air  at  14.70  psia  pressure  and  429.04JK,  373. 0°K,  and  288. 0°K 
temperature. 

b)  air  mixed  with  5.0%  steam  (by  weight)  at  14.70  psia  and  373. 0°K 
temperature.  The  results  are  presented  in  Figs.  3.1-3.13. 


3 . 2 Multistage  Compressor:  N.G.T.E.  #1 09 

The  design  details  for  this  compressor  can  be  found  in  Ref.  47. 

The  performance  of  the  compressor  was  calculated  for  operation  with 
the  following  working  fluids: 

(a)  air  at  14.70  psia  pressure  and  484.79°K 

(b)  air  mixed  with  10%  steam  (by  weight)  at  14.70  psia  pressure  and 
373. 0°K  temperature 

The  results  are  presented  in  Fig.  3.14. 
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3.3  P&W  TF-30  Compressor 


In  view  of  the  lack  of  consistent  design  and  experimental  data,  it 
has  not  been  possible  to  carry  out  performance  calculations,  although 
plans  have  been  made  to  carry  out  the  following  calculations: 


(a)  H.P.  compressor  performance  with  air  mixed  with  10%  by  v/eight 
steam. 

(b)  L.P.  compressor  performance  in  the  initial  stages  with  air 
mixed  with  10%  by  weight  water  in  uniform  droplet  form. 

(c)  Stage  10  with  sudden  evaporation  of  5%  water  and  with  air  mixed 
with  10%  by  weight  steam,  although  this  is  in  principle  equivalent 
to  the  calculation  described  under  3.1. 


4.  COMPRESSOR  AEROTHERMODYNAMICS 
The  compressor  aerothermodynamics  developed  here  is  based  on 
References  48-56.  However,  the  working  fluid  is  assumed  to  be  an  air- 
water  mixture  in  which  there  can  arise  a change  of  phase  from  liquid 
water  to  vapor.  In  view  of  the  fact  that  the  water  droplets  may  be 
subjected  to  gravitational  forces,  especially  in  inter-blade  spacings  and 
in  non-rotating  blade  passages,  gravitational  forces  ,.ave  been  included 
in  the  three-dimensional  two-phase  flow  equations. 

In  the  absence  of  gravitational  forces,  it  is  proposed  to  develop  a 
compressor  performance  calculation  method  based  on  the  axisynmetric  flow 
assumption. 


4.1  General  Considerations 
a)  State  Relations 

i)  GAS  PHASE:  The  gas  phase  is  considered  to  be  a mixture  of  air 
and  water  vapor  that  behaves  as  a perfect  gas.  Therefore, 


pg  = “g  Rs  \ 


(4.1) 


where 


Pg  = pressure  associated  with  the  gas  phase 
Rg  = gas  phase  gas  constant. 

ii)  LIQUID  PHASE:  The  liquid  phase  is  assumed  to  have  a constant 


density, 


pp  = 0.988  gm/crn3. 


(4.2) 


It  is  assumed  that  when  reaches  a value  of  0.40  for  tap  water  and 
air  (0.90  for  "pure"  water  and  air)  at  a pressure  of  one  atmosphere  a 
film  arises  on  the  casing  wall  (Section  2.5). 


b)  Stress-Strain  Relations 


l 


\ 


i)  Pressure:  The  mixture  pressure  (gas  continuum  + liquid 

"continuum")  is  assumed  to  be  equal  to  the  gas  phase  pressure, 
Pg.  It  is  assumed  that  the  water  droplets  do  not  contribute 
to  the  local  pressure  due  to  the  lack  of  randomness  in  local 
particle  motion. 

ii)  Viscosity:  A turbulent,  isotropic,  effective  coefficient  of 

viscosity  based  on  gas-phase  Reynolds  number  is  assumed 
throughout.  The  water  droplet  content  is  therefore  not 
accounted  for  in  estimating  the  mixture  viscosity, 

iii)  Viscous  force  on  droplets:  The  viscous  force  on  droplets  is 
based  on  a Reynolds  number  calculated  with  respect  to  the  gas 
phase  kinematic  viscosity  (effective),  the  relative  gas 
velocity,  and  the  local  statistical  mean  diameter  of  the  droplet 
(Section  2.6). 

iv)  Surface  tension:  The  droplet  breakup  is  based  upon  a critical 
value  of  Weber  number  (of  the  order  of  ten)  and  the  surface 
tension  for  water  is  taken  to  be 

a = 67.91  dynes/cm.  (4.3) 

c)  Water  Droplet  Fraction 

The  particulate  droplet  mass  fraction  is  defined  as  follows: 

mD 

°m  ' <yVolSi>/<VVoV  = ST  <4-4> 

r m 

where 

* 

am  = particulate  droplet  mass  fraction 
m = mass 
Vol  = Volume 

{ )-  = pertaining  to  the  mixture  as  a whole. 
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The  particulate  droplet  volume  fraction  is  defined  as  follows: 


ft 


jL 


av  = V01p/Vol-  . (4.5) 

Expanding  Equation  (4.4)  and  substituting  (4.5)  an  expression  is 
obtained  relacing  av  and  am 

°m  = % + pg>  (4'6) 

or 

o ~ a p /(p  - a [p  - p 1).  (4.7) 

v m Mg'  P mLMp  MgJ/  ' 

0^  is  a function  of  the  statistical  mean  diameter  and  number 

density  of  the  particulate  cloud. 


<r,  = | d 3 N 
v 6 s s 


(4.8) 


where  again 

dg  = the  Sauter  statistical  mean  diameter  of  a cloud  of 
droplets 

N$  = the  statistical  number  density,  i.e.  the  statistical 
number  of  ds  size  particles  per  unit  volume  of  overall 
mixture. 


4.2  Conservation  Equations 
a)  Assumptions 

i)  The  gas  phase  is  assumed  to  be  a perfect  gas  continuum, 
ii)  The  particulate  liquid  phase  is  assumed  to  be  a continuum  that 
does  not  contribute  to  the  local  pressure.  Two  exceptions  to 
this  continuum  assumption  are  the  determination  of  a)  the 
viscous  interaction  force  between  the  particulate  liquid  and 
gas  phases,  and  (b)  the  interfacial  heat  and  mass  transfers. 
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iii)  Mass  conservation  is  assumed  for  each  of  the  two  phases,  the 
gas  phase  and  the  particulate  liquid  phase,  allowing  for 
phase  change. 

iv)  The  liquid  mass  that  is  vaporir.ed  is  assumed  to  be  instantly 
and  intimately  mixed  with  the  gas  phase,  i.e.  large  scale 
vapor  concentration  gradients  and  the  momentum  of  the 
vaporized  mass  are  neglected. 

v)  On  the  scale  of  droplet  dimensions,  there  are  local  variations 
in  velocity,  temperature,  and  vapor  concentration.  The  evapora- 
tion of  individual  droplets  is  calculated  as  a process 
dependent  upon  such  local  properties, 

vi)  When  axial  synmetry  of  the  flow  is  assumed,  it  is  necessary  to 
treat  the  gravitational  body  force  as  being  negligible. 

| b)  Description  of  Coordinate  Systems 

The  equations  have  been  deduced  using  two  coordinate  systems, 
intrinsic  and  cylindrical, 
i)  Intrinsic  Coordinate  System  (Figure  4.1) 

The  m,e,n  intrinsic  system  is  a local,  right  hand,  orthogonal 
system.  The  m-di recti  on  lies  tangent  to  the  local  streamline 
direction,  and  the  e~direction  is  orthogonal  to  the  m-r  plane 
where  r is  the  radial  vector  defined  in  Figure  4.1.  The  n- 
direction  is  normal  to  the  m-di recti  on  in  the  local  meridional 
plane. 

ii)  Cylindrical  Coordinate  System  (Figure  4.2) 

The  cylindrical  system  is  fixed  to  the  axis  of  rotation  of  the 
compressor.  The  z-direction  lies  along  the  axis  of  rotation 
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in  the  direction  of  the  bulk  flow  and  r is  in  the  radial 


direction.  The  e-coniponent  is  in  the  circumferential  direc- 
tion defined  positive  with  respect  to  a right-hand  r,e, z 
orthogonal  coordinate  system, 
iii)  Rotation  of  the  Machine 

It  is  more  convenient  to  work  in  coordinate  systems  which  are 
stationary  with  respect  to  a rotating  blade  row.  Hence  fluid 
velocities  for  a rotating  machine  are  defined  as  follows: 

(1 ) Intrinsic  Coordinates 


V = W + U 


(4.9) 


where 


V = absolute  velocity  vector 
W = relative  velocity  vector 

0 = local  blade  rotational  velocity  vector  = + U0  + Un 

Um>0Q)0„  = the  m,6,n  components  of  0,  respectively, 
in  t?  n 


um " yu»*i»*2> 


(u ,<J>2 ) 


(4.10) 


Un  = fp|  (U  ><j>2  ) 


where 


U = the  magnitude  of  the  local  blade  rotational  velocity 
vector 

ti  ,Ua,U  = the  magnitudes  of  the  m,e,n  components  of  U 

iii  o n 

= the  angle  between  the  radius  vector  and  the  m-di  recti  on 
$2  = the  angle  between  the  radius  vector  and  the  n-di recti  on. 
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Therefore 


f 

= w„ 

+ 

Um 

V = 

w 

+ U 

gm 

gm 

m 

pm 

pin 

i 

= W A 

+ 

UA 

V = 

W A 

+ U 

ge 

ge 

e 

PO 

PO 

r 

= w 

+ 

U = 0 

V = 

w 

+ u 

gn 

gn 

n 

pn 

pn 

Again 


(2)  Cylindrical  Coordinates 


V = W + U, 


= 0 


where  similarly 


But 


Therefore 


u = ur  + Ue  + Uz. 


Ur  = Uz  = °* 


U = U 


e 


and  U = U.  = rw 

where 

to  = angular  velocity  of  the  rotating  blade  row. 
As  a result, 


(4.11) 


(4.12) 

(4.13) 

(4.14) 

(4.15) 


' 

= w 

V 

= w 

gr 

gr 

pr 

pr 

» 

= W „ + ru 

= w n 

ge 

ge 

pe 

pe 

f 

= w 

V 

= w , 

gz 

gz 

pz 

PZ 

+ rto  > 


(4.16) 


c)  Flow  Equations 

The  following  flow  equations  are  deduced  in  both  sets  of 
coordinate  systems  for  three-dimensional  as  well  as  axi symmetric 
flow.  They  consist  of  the  following: 
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i)  Mass  Conservation  equations 

ii)  Momentum  Conservation  equations 

iii)  Energy  Conservation  equations 

iv)  Variation  of  ay  with  respect  to  the  three  coordinate 
di rections 

v)  Radial  Equilibrium  equations 
The  equations  are  presented  in  Appendix  II  as  follows. 

11. 1.  Three-dimensional  Flow  Equations  in  Intrinsic  Coordinates 

11. 2.  Axisymmetric  Flow  Equations  in  Intrinsic  Coordinates 

11. 3.  Three-dimensional  Flow  Equations  in  Cylindrical  Coordinates 

1 1.4.  Axisymmetric  Flow  Equations  in  Cylindrical  Coordinates 


5.  DISCUSSION 


The  study  presented  in  this  Report  is  concerned  with  the  operation 
of  an  aircraft  compressor  with  an  air-water  mixture  ingested  into  the 
engine  either  from  puddles  on  runways  or  from  rain  at  altitude. 

The  water  in  the  air-water  mixture  may  be  found  in  several  states 
(discrete  droplets,  film  of  water,  evaporating  droplets  and  steam) 
between  the  inlet  and  outlet  of  a large  pressure  ratio  multi-stage 
machine,  such  as  the  TF-30  compressor.  Agglomeration  and  break-up  of 
droplets,  the  transport  processes  (mass  and  heat)  between  the  droplets 
and  the  surrounding  air  and  the  drag  of  droplets  are  not  phenomena  that 
are  well  understood  quantitatively.  They  are  central  to  gaining  an 
understanding  of  the  changes  in  the  state  of  water  in  the  system. 

The  presence  of  water  in  air  also  affects  the  compressibility  of 
the  mixture  and  therefore  the  local  Mach  number  at  any  location  in  the 
compressor  becomes  a function  of  the  state  and  content  of  water  in  the 
mixture. 

The  rotating  -flow  field  of  a compressor  has  a strong  influence  on 
the  state  and  distribution  of  water  since  centrifugal  forces  and  Coriolis 
forces  affect  both  discrete  droplets  and  films  of  water. 

Discrete  droplets  are  also  subject  to  gravitational  influence,  at 
least  at  low  rotating  speeds  and  in  inter- row  spacings  and  stators. 

In  examining  the  performance  of  a multi-stage  compressor  with  an 
air-water  mixture  flow,  it  seems  most  useful  to  divide  the  compressor 
into  groups  of  stages  on  the  basis  of  the  state  in  which  water  is 
expected  to  be  found  in  each  group.  In  doing  this,  there  arise  at  least 
two  important  uncertainties:  (1)  the  formation  of  a film  at  the  casing 

wall  and  (2)  the  gradual  evaporation  of  the  droplets. 
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Some  preliminary  investigations  conducted  on  the  performance  of  a 
single  stage  machine  with  air-water  droplet  mixture  flow  have  revealed 
that  the  water  droplets  cause  both  a reduction  in  mass  flow  as  well  as  a 
slight  change  in  efficiency,  mainly  due  to  incidence  and  Mach  number 


c? 


1/ ' 


Xl 


I?- 


effects . 

The  performance  of  stages  where  the  droplets  are  expected  to  be 
evaporating  cannot  be  estimated  with  any  degree  of  accuracy  based  on 
overall  stage  characteristics.  Nevertheless,  it  is  possible  to  argue 
that  the  stages  where  the  largest  amount  of  evaporation  occurs  must  show 
the  largest  change  in  performance  on  account  of  the  sudden  fall  in  air 
temperature  due  to  latent  heat  absorption  by  the  water. 

When  there  is  an  air-steam  mixture  flow,  it  is  possible  to  make  a 
rough  estimate  of  the  change  in  overall  stage  characteristics.  It  is 
found,  both  in  a single  stage  compressor  and  in  a multi-stage  machine, 
that  there  is  a noticeable  deterioration  in  performance  both  in  efficiency 
as  well  as  in  the  surgeline. 

In  order  to  perform  more  satisfactory  calculations,  it  is  imperative 
that  one  proceeds  to  set  up  the  equations  of  aerothermodynamics  for  air- 
water  mixture  flow  in  a compressor.  This  task  has  been  completed  both  in 
general  three-dimensional  coordinates  as  well  as  in  axisymmetric 
coordinates.  Such  equations  involve  several  quantities  that  can  only  be 
estimated  very  roughly  and  hence  can  perhaps  be  treated  only  in  a para- 
metric form:  (1)  volume  fraction  of  the  water  droplets,  (2)  statistical 

mean  diameter  and  number  density  of  droplets,  (3)  mass,  momentum  and 
energy  transfer  between  the  droplets  and  the  air,  and  (4)  degree  of  mixing 
between  evaporated  water  and  the  surrounding  air. 
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The  equations  for  air-water  mixture,  when  compared  with  the  standard 


compressor  flow  equations  for  air  flow,  reveal  several  ways  in  which 
computational  programs  set  up  for  air  flow  can  be  adapted  for  use  with 
air-water  mixture  flow. 

In  examining  the  problem  from  the  point  of  view  of  alleviating  the 
effects  of  v/ater  ingestion,  one  can  arrive  at  several  methods  of  attack 
on  the  problem: 

(a)  Removal  of  as  much  water  as  possible  at  the  casing  walls. 

(b)  Bleed-in  of  air  into  the  latter  stages  of  the  compressor  from 
another  stream. 

(c)  Variable  geometry  engine  with  adjustable  air  flow  paths. 


5.1 . Reconroendations 

The  recommendations  may  be  divided  broadly  into  two  categories: 

(i)  Further  diagnostic  studies. 

(ii)  Future  studies. 


5.1.1.  Further  diagnostic  studies 

(1)  Estimation  of  performance  of  selected  multi-stage  machines 
(incorporating  bleed-in  and  bleed-out)  on  the  basis  of  stage  character- 
istics when  there  is  a simultaneous  change  in  temperature  from  the  design 
value  and  an  addition  of  steam. 

(2)  Estimation  of  the  effect  of  bypass  ratio  on  the  performance  of 
compressors  with  two-phase  flow. 

(3)  Study  of  two-phase  flow  in  a diffusing  passage. 

(4)  Adaptation  of  simple,  existing  computational  programs  for  air 
flow  to  air-water  mixture  flow  on  the  basis  of  parameterized  droplet 
dynamics  and  transport  processes  for  the  study  of  a selected  multi-stage 
machine. 


(5)  A flow  visualization  study  on  a multi-stage  compressor  to 
establish  the  water  droplet  dynamics  when  water  is  ingested  into  the  air 
stream  at  entry  to  the  compressor 

(6)  A flow  visualization  study  on  a multi-stage  compressor  with 
selected  inlet  configurations  to  establish  the  water  droplet  dynamics 
when  water  is  injected  into  the  air  stream  at  entry  to  the  inlet. 

(7)  A flew  visualization  study  on  a rotating  tire  in  a wind  tunnel 
equipped  with  a moving  belt  floor  to  establish  the  influence  of  a few 
selected  tire  parameters  on  the  jet  created  by  the  motion  of  the  tire. 

5.1.2.  Future  studies 

(1)  Development  of  a comprehensive  computational  program  for  air-water 
mixture  flow  through  an  axial  compressor  with  appropriate  droplet  dynamics 
and  transport  coefficients. 

(2)  Analysis  of  the  performance  of  compressors  using  the  computational 
program  developed  to  establish  the  scope  and  implications  of  variable 
cycle/geometry  schemes. 

(3)  (a ) Study  of  bleed-out  of  water  in  early  stages. 

(b)  Study  of  bleed-in  of  air  into  the  latter  stages. 

(4)  Determination  of  classes  of  compressors  to  which  different 
classes  of  controls  can  be  applied. 

(5)  Study  of  the  effects  of  engine-airframe  integration  parameters 
on  the  utilization  of  different  kinds  of  controls. 

j (6)  Study  of  the  water  jet  created  by  a moving  aircraft  tire  on  a 

1 

1 

J wet  runway  in  relation  to  the  location  of  the  engine  inlet. 

1 - 

(7)  Determination  of  the  effect  of  water  vapor  on  after-burner 
stability,  particularly  in  relation  to  the  fuel  vaporization  and  combus- 
tion characteristics. 
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FIGURE  3.31 
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Three-Dimensional  Intrinsic 
Coordinate  System 


Axisymmetric  Intrinsic 
Coordinate  System 
Wennerstrom 
(Ref.  56) 


FIGURE  4.1 

INTRINSIC  COORDINATE  SYSTEM 
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FIGURE  4.2 
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Estimate  of  Water  Lifted  by  Wheels 
An  estimate  of  the  amount  of  water  thrown  up  by  an  aircraft  wheel 
from  a puddle  of  depth  equal  to  h can  be  obtained  as  follows. 

Let  the  direction  of  motion  of  the  wheel  be  x and  the  transverse 
direction  on  the  ground  be  y.  Refer  to  Figure  A.I.l. 

It  has  been  pointed  out  (Ref.  57)  that  the  relative  flow  between 
the  wheel  and  the  water  is  hypersonic  in  the  sense  that  the  Froude  number, 
U /(g  h)z,  is  of  the  order  of  100  when  h is  of  the  order  of  a cm.  in 
depth,  and  where  Uc  is  the  speed  of  the  undercarriage  of  the  wing  and 
hence  of  the  order  of  tens  of  meters  per  second.  Accordingly,  one  can 
postulate  (on  the  basis  of  blast  wave  analogy)  a bow  wave  or  hydraulic 
jump,  as  shown  in  F*g.  A.I.l,  at  a distance  yg  from  the  wheel. 

Now,  a jet  of  spray  can  be  assumed  from  the  wheel  (Ref.  58)  on 
account  of  the  action  of  the  relative  wind  on  the  surface  of  the  moving 
water.  The  time  rate  of  "pumping"  of  the  v/ater  can  be  written  as 


where  p is  the  density  of  air,  p the  density  of  v/ater,  yR  the  time 
derivative  of  y^,  and  c*  a constant  equal  to  approximately  one. 

In  order  to  obtain  yg,  consider  the  momentum  of  the  mass,  pw  yg  h, 
contained  in  a strip  of  v/ater  of  unit  width  and  the  frictional  resistance 
offered  by  the  ground,  fp  yg  yg2,  where  f is  the  friction  coefficient, 
of  the  order  of  10  . The  two  forces  must  balance  each  other  together 
with  the  momentum  loss  in  the  spray.  It  can  then  be  shewn  that 


1 n V2 

- tf1  - r> 


(A. 1.2) 
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where 


and 


n,  - 0.5  Ci(Pa/Pw) 


H 


+ f 
ni 


(A. 1.3) 
(A.I.4) 


Substituting  Eg.  (A. 1.1)  into  Et.  (A. 1.2),  the  amount  of  water 
sprayed  behind  the  wheel  per  unit  distance  of  wheel  travel  can  be  written 
as  follows. 

- VW*  h2/ci-  (fl-I-5> 

Thus,  the  mass  of  water  thrown  up  from  a puddle  0.5  meters  long  and  2 cm. 
deep  is  of  the  order  of  10  kg.  during  a conventional  take-off.  This  does’ 
not  yet  account  for  the  water  sprayed  forward  of  the  wheel  in  its  jet 
stream.  During  a rain  storm,  over  a runway  length  of  1,000  meters  with  a 
1 cm.  depth  of  water  on  the  runway,  the  amount  of  water  lifted  off  the 
surface  can  be  as  much  as  several  metric  tons.  Motion  pictures  of  wheel 
sprays  have  shown  rather  large  pumping  characteristics. 

An  estimate  of  the  spray  envelope  from  the  wheel,  y^  and  zB  (the 
height  reached  by  the  spray),  is  rather  more  difficult  to  estimate.  Some 
experimental  studies  on  a typical  wheel  in  a moving-belt  equipped  wind 
tunnel  is  indicated  before  proceeding  too  far  in  the  analysis. 
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b;  Momentum  "Miervation  Equations: 
Gas  Phase 
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c)  Energy  Conservation  Equations: 
Gas  Phase 
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d)  Variation  of  cTv  with  respect  to  the  three  coordinate  directions: 
m-di recti  on 

The  equation  is  obtained  as  follows.  The  e-component  gas  phase 
momentum  conservation  equation  is  differentiated  with  respect  to  e.  Next 
the  m-component  gas  phase  momentum  conservation  equation  is  solved  in 
terms  of  ( X - crv).  These  two  equations  are  then  combined  with  the  gas 
phase  mass  conservation  equation  to  obtain  the  final  result. 


(A.II.ll) 

e-direction 

The  m-component  gas  phase  momentum  conservation  equation  is  differ- 
entiated with  respect  to  m.  Next  the  e-component  gas  phase  momentum 
conservation  equation  is  solved  in  terms  of  ( X - CTV  ).  These  two 
equations  are  then  combined  with  the  gas  phase  mass  conservation  equation 
to  obtain  the  final  result. 
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The  equation  is  obtained  by  first  differentiating  the  n-component 
gas  phase  momentum  conservation  equation  with  respect  to  n.  Then  the 
o-component  gas  phase  momentum  conservation  equation  is  solved  in  terms 
of  ( 1 - <rv  ) and  combined  with  the  differentiated  n-component  gas 
phase  momentum  conservation  equation  to  obtain  the  final  result. 
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e)  Equilibrium  Equations: 

Gas  Phase 

An  enthalpy-entropy  relation  valid  for  the  gas-phase  flow  can  be 
written  as  follows. 
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(A. II. 14) 


where 


h 


5 


= enthalpy  of  the  gas  phase  per  unit  mass  of  the  gas  phase 
(joules/Kgg) 

= rise  in  entropy  of  the  gas  phase  (per  unit  volume  of  mixture) 
due  to  the  action  of  body  and  surface  forces  on  the  gas  phase, 
excluding  those  forces  imposed  by  the  liquid  phase  (joules/°K 
Volj) 
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= mass  averaged  gas  phase-liquid  phase  temperature  (°K) 

Js-j  = rise  in  entropy  of  the  gas  phase  (per  unit  volume  of  mixture) 
due  to  the  action  of  forces  on  the  gas  phase  imposed  by  the 


liquid  phase  (joules/°K  Vol - ) 


= the  heat  loss  from  the  gas  phase  (per  unit  volume  of  mixture) 
to  the  surroundings,  excluding  the  heat  lost  directly  to  the 


liquid  phase  (joules/Vol-) 


dc^  = heat  loss  from  the  gas  phase  (per  unit  volume  of  mixture) 


directly  to  the  liquid  phase  (joules/Vol-)' 


By  definition 


V»v  V3 


(A. 11.15) 


where 


}-jr  = total  enthalpy  of  the  gas  phase  per  unit  mass  of  the  gas  phase 

o 

(joules/Kgg). 

Combining  (A. 11.14)  and  (A. II. 15), 

d ^ d Hj  “ d - V^©  d V^o  -*  1 x 

* J (A. II. 16) 

- Tj  -i “ Ty  d + dj. 

Combining  (A. 11.16)  with  each  of  the  component  gas  phase  momentum 
equations  (eliminating  the  pressure  gradient  term)  yields  the  following 
equilibrium  results. 
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(A. 11.19) 


Liquid  Phase 


An  internal  energy-entropy  relation  valid  for  the  liquid  phase  flow 
can  be  written  as  follows. 

dufs  JL  Tp  d5p  -7cjp  cU*  “ cjc^  n cl  1 (A. 11.20) 

crv  P/j  J 


where 


Up  = internal  energy  of  the  liquid  phase  per  unit  mass  of  the  liquid 
phase  (joules/Kgp) 

dsp  = rise  in  entropy  of  the  liquid  phase  (per  unit  volume  of  mixture) 

due  to  the  action  of  body  and  surface  forces  on  the  liquid  phase, 

excluding  those  forces  imposed  by  the  gas  phase  (joules/°K 

Vol-) 

m 

cjc^p  = heat  Toss  from  the  liquid  phase  (per  unit  volume  of  mixture) 
to  the  surroundings,  excluding  the  heat  lost  directly  to  the 
gas  phase  (joules/Vol-). 
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By  definition 


Up  ~ Up  -r  Vp ^ 


(A. 11.21) 


where 


U = total  internal  energy  of  the  liquid  phase  per  unit  mass  of  the 
v P 

liquid  phase  {joules/Kgp) . 

Combining  (A. 11.20)  with  (A. 11.21)  and  expanding 
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(A. 11.22) 


"7"p  d 5 p “ "Tfj  p cU  J - cl  p + Q Q|^ 


Combining  (A. 11.22)  with  each  of  the  component  liquid  phase  momentum 
conservation  equations  yields  the  following  equilibrium  results. 
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II. 2.  .,^'mmt.ric  flow  EquationsJnJntHnsic  Coordinates 
a)  Mass  Conservation  Equations: 

Gas  Phase 


if  JL  trLl'^^WynV  Qhz 

r 


= hr>x 


Liquid  Phase 


1 U(r<r1,ppWPJ+  r<r^p  iV^ 


r Ld  vn 


1 V-*  4 


m y 


88 


iL 


(A. II. 25) 


(A. 11.26) 


(A. II. 27) 
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b)  Momentum  Conservation  Equations: 
Gas  Phase 
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c)  Energy  Conservation  Equations: 
Gas  Phase 
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Liquid  Phase 
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d)  Variation  of  crv  with  respect  to  the  three  coordinate  directions: 
ro-direction 


6-direction 
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e)  Radial  Equilibrium  Equation: 

These  equations  are  deduced  utilizing  the  method  of  Wennerstrom 
(Ref.  56). 
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Liquid  Phase 
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It  may  be  pointed  out  here  that  in  general  three-dimensional  flow, 
it  is  necessary  to  set  up  equilibrium  conditions  in  all  of  the  three 
directions  (Eqs.  A. II. 17,  A. 11.18,  A. II. 19  and  A. 11.23,  A. 11.25,  A. 11.25). 

I1,3*  lilTse-dimenslonal  Flow  Equations  in  Cylindrical  Coordinates 
a)  Mass  Conservation  Equations: 

Gas  Phase 
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Liquid  Phase 
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b)  Momentum  Conservation  Equations  : 
Gas  Phase 
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Liquid  Phase 
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d)  Variation  of  CTv  with  respect  to  the  three  coordinate  directions: 
r-direction 

The  equation  is  obtained  as  follows.  The  e-component  and  the  z- 
component  gas  phase  momentum  conservation  equations  are  differentiated 
with  respect  to  8 and  z,  respectively.  Next,  the  r-component  gas  phase 
momentum  conservation  equation  is  solved  in  terms  of  ( 1 -c rY  ).  These 
three  equations  are  then  combined  with  the  gas  phase  mass  conservation 
equation  to  obtain  the  final  result. 
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e-direction 

This  equation  is  obtained  as  follows.  The  r-component  and  the  z- 
component  gas  phase  momentum  conservation  equations  are  differentiated 
with  respect  to  r-  and  z,  respectively.  Next,  the  e-component  gas  phase 
momentum  equation  is  solved  in  terms  of  ( l-o~v  ).  These  three  equations 
are  then  combined  with  the  gas  phase  mass  conservation  equation  to  obtain 
the  final  result. 
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This  equation  is  obtained  as  follows.  The  r-component  and  the  e- 
compor  nt  gas  phase  momentum  conservation  equations  are  differentiated 
with  respect  to  r and  0,  respectively.  Next  the  z-component  g -s  phase 
momentum  conservation  equation  is  solved  in  terms  of  ( \ ~JTV  ).  These 
three  equations  are  then  combined  with  the  gas  phase  mass  conservation 
equation  to  obtain  the  final  result. 


+ Wyr  JL 

P-jV/jX  Jr 


I li  .y  t I 


Vy  if  i,  - \jz  t V-  a i^\  S iy^V 

„ Jr  \ > ) jo 


- 1 


-i-kj ajxtj  -r  Wei  iki-]  + C r~«yL  -t  d-.j 

•'  -»*  J L *J 

V J 


1..  \ 0 sUi  | 

j r I r / J j l J 

J Pj  i Vvi£_)  ~t  'j-  (0  V’o  -+•  -i.  ( ) O..  \v..i  ) IT  pc  War  5 

L o v J o J i J »• 


(A. 11.53) 


War  Q Wr 
Jr 


+ * u‘ ) ^-Vy _ 

\ r J J-i> 


J Wyi  1 1 T £ Wy  1 X 


Wr.r 

' j r 


100 


r 


I 

1 


+ ( Wfrg  -t  S Vat  - Wftdtru^)  -<■  V/f.i  ^U4r  1 " Pc  ( WcSiri-')  x 

\ ^ ilt*  r It  1 n 1 

| *f  Wy  ( | -4  W}2.  S Wf|6  1 

\ r J ^ d 1 r / ^2  J 


JL 

de 


Warf  o W^a  -+  U)  ^ *t  f 

\ Jr  J 


-2 

V/^r  J SWr  r J 1 f j ^ ^Ag  -t  Wjr/  W»-?  + u>  j •V\/v/c' 7.  SV/^a  I 

l d'  / l r /a®  \ r J J ^ J 


e)  Equilibrium  Equations: 

Gas  Phase 

These  equations  are  obtained  by  combining  (A. I I. 16)  with  each  of  the 
three  component  gas  phase  momentum  conservation  equations. 
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Liquid  Phase 

These  equations  are  obtained  by  combining  {A. II. 22)  with  each  of  the 
three  component  liquid  phase  momentum  conservation  equations. 
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n-4-  Axi symmetric  Flow  Equations  in  Cylindrical  Coordinates 
a)  Mass  Conservation  Equations: 
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b)  Momentum  Conservation  Equations: 


Gas  Phase 


r-component 


( j Qr,  r Wc  v $ \\A  «"  — ( \aA  J ~t  |*U  ) ■+-  Wq  7 O lv  r, 

‘ ()  •)  \ W J j . ' 


r 

— o 


iA  ■*  Fn.v  + A'.,  ■*  ;ri,- 

\ j 4 

dr 


e-component 


( J ~ 3~y  VA  ^ S Vv^c>  ■+  lo  ^ t W<^  r | Wc_3  * l- 


+ H*  S Wc^e  ] - F^y  r-./;  e>  t 


z-component 


( 3.  ~<T', y.  ^ \\/c^r  + VlAj  2.  S -2,  — 

c)r  . 


"*■  Fa*  i ^ I + 


Liquid  Phase 


r-component 


-Uairv)  + Vp  z *5  W k>  r 


iL 


~ /: 


Bpr  t )'Y^  Fj t 


(A. 11.65) 


e -component 


c r„ 


fp 


^Pr  ( S Wpa  -tU> 
c)r 


■*  Vv^jOir  | Wpa  + u> 

\ y*  / 


4 SlWs? 

IT 


o>  * ^ ^pe  ~ ^Xe 


(A. II. 66) 


^component 


0~\ 


Hr 

c)  2 


(A. 11.67) 


c)  Energy  Conservation  Equations; 
Gas  Phase 
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Liquid  Phase 
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d)  Variation  of  <rv  with  respect  to  the  three  coordinate  directions: 
r-di recti  on 

- - lUv  ■+  S_  1 -t  It-z.  " 

<d  r Wy  Wy  d 2 L d i 


n 


XJy  ft  z -*  V/y  S Wy 

^■V"  <}  2.  J 


i ^ Bc.  y Fvn  ) -r  . ft) 

1 \ I i J 


dr. 


Vy  SVy  •+  Wy  S r — 1 W-^ dif  lO ) 

c)  r <)  i r 


-1 


5L  { Vy  ) -t  Sir  \v(.  z.) 
dr  j z 


rp<  Y/«*  JL 


'ft  o 


dz. 


(A. 11.70) 

Vy  ^ Vv/Oj  <L  *ft  Wy 

dr  ^ z. 


-1 


Wy  u W«j  z -v  Vi/pj  z S Wy 
dr*  d 2 J 


f 

) 
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9-direction 


So  y/  - O 

J © 


z-di recti  on 


S<rv  = S~  i > 

J-  z p^z.  fjH2-  ^ 


2.  J 

V\/or  *6  Wo.  V ~ ( \v4qH  l1-'  ) -i  We.  £ *S>  W-\  1 1 -}•  1 — 

L 3 T ~V-  'ati  rft'Wy 


fA  -,  /v  1 FW  - Sp*  Vnr  &.t  ->  WW  ^J4-z  * 

^ z J L dr  C/ZJ 


p^  £.(r\V«jr)  -»  X 
Jr  c)2 


yp<  W c r 

! <)  >5 


r z . 

£ W^v  ^ Wyr  - l W<^  o ~t  1 1/* ) -*  vv/tj  i y 

J r [ <W  r J 2 - 


B-component 


,pf  f Wr  * V/f ^ * 2 w 


- f a ^ 'I'Vae  " ^1 


z- component 


CT  \f  |Pp  V p f 


-<bW1o< 
^ Z 


-tvf  a-+ru>  ) S Wj3d 

^ Z 


= - (Tv  Pf  12*  + ^ W -V  ^ "4* 

lr^z  (Ji  a d7- 


<)  z 


if  U,  ',.  Vr.CfWIlf  PRIKIINO  OtflCf:  !9Jfc  — ?W-501/» 
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